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Abstract. We investigate the thermal and kinematic configuration of a sunspot penumbra using very high spectral and spatial 
resolution intensity profiles of the non-magnetic Fei 557.6 nm line. The dataset was acquired with the 2D solar spectrometer 
TESOS. The profiles are inverted using a one-component model atmosphere with gradients of the physical quantities. From 
this inversion we obtain the stratification with depth of temperature, line-of-sight velocity, and microturbulence across the 
penumbra. Our results suggest that the physical mechanism(s) responsible for the penumbral filaments operate preferentially 
in the lower photosphere. We confirm the existence of a thermal asymmetry between the center and limb-side penumbra, the 
former being hotter by 100-150 K on average. We also investigate the nature of the bright ring that appears in the inner penumbra 
when sunspots are observed in the wing of spectral lines. It is suggested that the bright ring does not reflect a temperature 
enhancement in the mid photospheric layers. The line-of-sight velocities retrieved from the inversion are used to determine the 
flow speed and flow angle at different heights in the photosphere. Both the flow speed and flow angle increase with optical 
depth and radial distance. Downflows are detected in the mid and outer penumbra, but only in deep layers (logTsoo < —1.4). 
We demonstrate that the velocity stratifications retrieved from the inversion are consistent with the idea of penumbral flux tubes 
channeling the Evershed flow. Finally, we show that larger Evershed flows are associated with brighter continuum intensities 
in the inner center-side penumbra. Dark structures, however, are also associated with significant Evershed flows. This leads us 
to suggest that the bright and dark filaments seen at O'.'S resolution are not individual flow channels, but a collection of them. 
Our analysis highlights the importance of very high spatial resolution spectroscopic and spectropolarimetric measurements for 
a better understanding of sunspot penumbrae. 
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1. Introduction 

The fine structure of the penumbra is intimately connected 
to the Evershed flow, the most conspicuous dynamical phe- 
nomenon observed in sunspots (e.g., Solanki 2003; Thomas & 
Weiss 2004). Therefore, a proper observational characteriza- 
tion of this fine structure is essential to understand the Evershed 
flow and, more generally, the nature of the penumbra itself. 

Studying the fine structure of the penumbra is difficult be- 
cause of the very small scales involved. Recent imaging ob- 
servations with the Swedish 1-m Solar Telescope (Scharmer 
et al. 2003; Rouppe van der Voort et al. 2004) have revealed 
that even at a resolution of ff.' 1 some of the filaments that form 
the penumbra may be unresolved. That is, different structures 
(each having different properties) may coexist in the resolu- 
tion element, and we are not able to separate them. A proper 

* Operated by the Association of Universities for Research in 
Astronomy, Inc. (AURA), for the National Science Foundation 



characterization of the penumbra calls not only for diffraction- 
limited imaging observations, but also for more quantitative 
spectroscopic and spectropolarimetric analyses at the highest 
resolution possible. By interpreting the shape of spectral lines 
we can derive the thermal, magnetic, and kinematic configu- 
ration of the penumbra via the Zeeman and Doppler effects. 
Unfortunately, current spectroscopic observations do not attain 
resolutions better than ~ O'.'S, and this limit is often reached at 
the cost of poor spectral resolution. The situation is even worse 
in the case of spectropolarimetric measurements: the best an- 
gular resolution of existing solar polarimeters is about 1". 

This paper is the third of a series devoted to the anal- 
ysis of 2D spectroscopic observations of a sunspot. The 
data were acquired with the TElecentric SOlar Spectrometer 
(TESOS; Kentischer et al. 1998; Tritschler et al. 2002) and 
the Kiepenheuer Adaptive Optics System (KAOS; Soltau et 
al. 2002; von der Liihe et al. 2003) at the German Vacuum 
Tower Telescope of Observatorio del Teide (Tenerife, Spain). 
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These observations combine high spectral (A/AA ~ 250000) 
and spatial (075) resolution (Tritschler et al. 2004; hereafter re- 
ferred to as Paper I). In Paper II of this series (Schlichenmaier 
et al. 2004), we investigated the properties of the Evershed flow 
through a bisector analysis of the observed intensity profiles. 

Here we perform a full inversion of the same data set in 
order to study the thermal and kinematic configuration of the 
penumbra. Our analysis allows us to explain a number of find- 
ings reported in Paper I and to remove some of the uncertainties 
associated with the methods employed in Paper II. At the same 
time, we gather new information on the thermal properties of 
the penumbra. This is particularly important because very few 
thermal studies of sunspot penumbrae have been published 
to date (del Toro Iniesta et al. 1994; Balasubramaniam 2002; 
Rouppe van der Voort 2002; Sanchez Cuberes et al. 2005). 

A brief account of the observations and details of the inver- 
sion procedure are given in Sect. [2 In Sect. [3] we discuss the 
maps of physical quantities inferred from the inversion. The 
thermal and kinematic configuration of the penumbra is de- 
scribed in detail in Sects. |4] and [5] Some implications of the 
results are presented in Sect.|6j and a summary of our findings 
is given in Sect.0 

2. Observations and data analysis 

2.1. Observations 

On July 5, 2002, TESOS was used to measure the intensity pro- 
files of Fe i 557.6 nm in the main sunspot of NOAA AR 10019 
(cf. Paper I). At the time of the observations, the spot was lo- 
cated 23° off the disk center. The line was scanned at 100 wave- 
length points with a wavelength step of 8.4 pm, corresponding 
to overcritical spectral sampling. The time needed to complete 
the scan was 37 s. Operated in high resolution mode, TESOS 
yielded a pixel size of 07089 x 07089. During the scan, the 
KAOS system provided real-time correction of wavefront dis- 
tortions due to turbulence in the Earth's atmosphere. 

2.2. Diagnostic capabilities of Fe i 557.6 nm 

We use the intensity profiles of Fei 557.6 nm to derive the 
stratification of temperature and velocity in the penumbra. 
Therefore, it is important to examine the properties of this line 
in some detail. Fei 557.6 nm is a fairly strong photospheric 
line with zero effective Lande factor, i.e., it does not undergo 
any Zeeman splitting. The line is often quoted to be tempera- 
ture insensitive, but this is too a simplistic statement. Figure 5 
of Paper I shows that the equivalent width of Fei 557.6 nm 
changes by non-negligible amounts when going from the quiet 
sun to the umbra, where the temperatures are rather different. 

In the following, we elaborate on the temperature sensi- 
tivity of Fe i 557.6 nm by means of response functions (RFs; 
see Ruiz Cobo & del Toro Iniesta 1994 and del Toro Iniesta 
2003). The upper panel of Fig.^shows the RF of Fe 1 557.6 nm 
to temperature perturbations as evaluated in the penumbral 
model atmosphere of del Toro Iniesta et al. (1994; hereafter 
TTR). This plot tells us how much the specific intensity at A 
changes when the temperature of the model is increased by 
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Fig. 1. Response functions of the Fe i 557.6 nm line to tem- 
perature (top) and line-of-sight velocity (bottom). The two 
plots indicate the percentage by which the specific intensity 
at A changes when the temperature (resp. LOS velocity) is in- 
creased by 1 K (resp. 1 m s -1 ) in a layer of optical width 
AlogT5oo =0.1 and optical depth Tsoo- The thermal stratifi- 
cation of the penumbral model of TTR has been employed for 
these calculations. 



1 K at a given optical depth t$oo- The figure demonstrates that 
the line reacts to temperature changes. The sensitivity of any 
spectral line to temperature perturbations is determined by two 
competing effects: changes in the absorption and in the emis- 
sion. Usually, only opacity variations are considered (e.g., Gray 
1988, Chapter 13) but, as pointed out by Cabrera Solana et al. 
(2005), the dominant contribution is that of the source function. 
Neglect of this dominant factor is the likely cause of the line 
being quoted as temperature insensitive in most publications. 

Figure |2 shows the temperature dependence of the equiv- 
alent width for a set of neutral iron lines commonly used in 
solar physics (adapted from Cabrera Solana et al. 2005). It is 
apparent from this figure that Fe i 557.6 nm is one of the most 
temperature sensitive lines. Indeed, it shows a larger sensitivity 
than lines presumed to be appropriate for temperature measure- 
ments such as Fei 524.7 and 525.0 nm. The large variation of 
the equivalent width of Fe i 557.6 nm with temperature is con- 
sistent with the observed behavior of the line in the quiet sun, 
umbra and penumbra (Fig. 5 of Paper I). 

Let us now turn our attention to the sensitivity of Fe 1 557.6 
nm to line-of-sight (LOS) velocities. The bottom panel of Fig.^ 
shows the RF of the intensity to velocity perturbations as eval- 
uated in the penumbral model of TTR. These RFs correspond 
to a model where all velocities have been set to zero, so they 
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Fig. 2. Variation of the equivalent width W of 9 neutral iron 
lines when the temperature is increased by 1 K in all photo- 
spheric layers. The calculations have been carried out using the 
thermal stratification of the penumbral model of TTR. 

are symmetric about the line center (models in which the LOS 
velocity varies with height would result in asymmetrical RFs 
and slightly different sensitivities in the different atmospheric 
layers). As can be seen, the line is sensitive to velocities in a 
height range from about logrsoo = to -3.5. The maximum 
sensitivity is attained in the mid photosphere. In the deep and 
high layers where the continuum and the line core are formed, 
respectively, the sensitivity to velocities is much reduced. The 
reason has been explained by Cabrera Solana et al. (2005): at 
a given wavelength, the intensity variation induced by Doppler 
shifts is proportional to the slope of the line at that position. 
Thus, the intensity changes are maximum near the line wings 
(which sample the mid photosphere) and minimum near the 
continuum and the line core, where the slope of the intensity 
profile (dl/dA) quickly goes to zero. The lower panel of Fig.^ 
also demonstrates that the intensity observed at a fixed wave- 
length is sensitive to LOS velocities in a broad range of optical 
depths. That is, the Doppler shifts measured from line bisectors 
at a given intensity level contain information from a significant 
fraction of the atmosphere. Thus, it is difficult to ascribe them 
to any particular layer. This limitation of bisector analyses has 
long been recognized (Maltby 1964; Rimmele 1995). 

2.3. Inversion procedure 

The intensity profiles observed with TESOS have been inverted 
in order to determine the stratification of temperature, LOS ve- 
locity, and microturbulence with optical depth. To this end, 
we have used the SIR code (Ruiz Cobo & del Toro Iniesta 
1992). Despite the high spectral resolution of our observations, 
the number of data points does not allow us to use complex 
atmospheric models. Hence, we have adopted a simple one- 
component model with gradients of the physical parameters. 

We derive the temperature stratification by modifying the 
thermal structure of the penumbral model of TTR with two 
nodes. For a better recovery of this parameter, the observed 
profiles are normalized to the continuum intensity of the av- 
erage quiet sun profile. Such an absolute normalization implies 
that the continuum of the individual profiles carries information 



about the temperature in the deep atmospheric layers. The two 
nodes used to infer the thermal structure allow for changes in 
the gradient and absolute position [e.g., T(tsoq = 1)] of the ini- 
tial temperature stratification. We assume local thermodynami- 
cal equilibrium (LTE). Since non-LTE effects may be important 
in the upper layers, our temperatures are reliable only up to, 
say, logT5oo ~ -3. Electron pressures are computed from the 
inferred temperatures using the equation of hydrostatic equilib- 
rium and the ideal gas law. 

The inversion also yields the stratification with depth of the 
LOS velocity and microturbulent velocity. In both cases we as- 
sume linear variations with logrsoo, i.e., two nodes for each 
parameter. In the case of the LOS velocity, this choice seems 
to be appropriate considering that most observed bisectors pos- 
sess quite linear shapes (Paper II). From the inversion we also 
determine the (height-independent) macroturbulence needed to 
reproduce the observed line widths. The synthetic profiles are 
convolved with the instrumental profile of TESOS before being 
compared with the observed ones. To describe the broadening 
of the line due to collisions with neutral hydrogen atoms we 
use the quantum mechanical formulation of Anstee, Barklem, 
& O'Mara (see Barklem et al. 2000). The total number of free 
parameters is 7, which compares favorably with the number of 
observables (100). 

Without additional spectral lines observed, or measure- 
ments of the four Stokes parameters, it is not possible to esti- 
mate the amount of stray light contamination. We expect, how- 
ever, that stray light arising from seeing fluctuations is much 
reduced with respect to earlier observations that did not benefit 
from adaptive optics. By contrast, parasitic light inside TESOS 
has been corrected for. In Paper I we determined that ~ 5% of 
the mean continuum intensity contributes to each filtergram in 
the form of ghosts and scattered light. This parasitic light has 
been subtracted from the individual images before extracting 
the intensity profiles. 

Our inversion approach avoids important shortcomings of 
simpler analyses. Since the complete radiative transfer prob- 
lem is solved, we effectively separate the contribution of tem- 
peratures and LOS velocities to the observed profiles using the 
information provided by the RFs. From Fig. \l\ it is clear that 
the intensity at a fixed wavelength depends on the temperature 
and LOS velocity stratification along the atmosphere. Thus, 
we determine both parameters simultaneously in order not to 
misinterpret thermal effects as velocity effects and vice versa. 
Another important improvement is that we can ascribe the in- 
ferred velocities and temperatures to specific optical depths, be- 
cause our procedure takes into account the finite width of the 
intensity contribution functions. 

The simple one-component model atmosphere used in this 
investigation has several limitations. First, it does not allow for 
any unresolved structure in the pixel. Although our spatial res- 
olution is one of the highest ever reached in spectroscopic stud- 
ies, we are still far from resolving the individual constituents of 
the penumbra. Spectropolarimetric measurements suggest that 
at least two magnetic components are necessary to understand 
the penumbra (Schlichenmaier & Collados 2002; Bellot Rubio 
et al. 2003, 2004; Borrero et al. 2004, 2005). One of the two 
components carries most of the Evershed flow, while the other 
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Fig. 3. Examples of observed (dots) and best-fit (solid lines) in- 
tensity profiles emerging from the center-side (left) and limb- 
side (right) penumbra. The vertical dashed lines indicate the 
position of zero velocities (see Paper I). Note the strong asym- 
metries of the line, with an extended blue wing in the center- 
side penumbra and an extended red wing in the limb-side 
penumbra. 

is essentially at rest. If the flow channels are horizontally unre- 
solved and we do not account for this possibility, the LOS ve- 
locities determined from the inversion would represent lower 
limits to the real velocities. Second, we have chosen a very 
simple functional dependence to describe the run of the LOS 
velocity with height. In an uncombed penumbra consisting of 
penumbral flux tubes embedded in a background magnetic at- 
mosphere, lines of sight that pierce the tubes would 'see' dis- 
continuities in the LOS velocity as they encounter the upper 
and lower boundaries of the tubes. Our assumption of linear ve- 
locity stratifications would provide only a very rough descrip- 
tion of such discontinuous velocity stratifications. In that case, 
it would be the magnitude and sign of the LOS velocity gra- 
dient inferred from the inversion which would inform us about 
the velocity inside the flux tubes and the height position of the 
flow channels. 

Despite these shortcomings, the adopted model atmosphere 
does an excellent job in explaining the observed intensity pro- 
files. On average, the residuals of the fit are only slightly larger 
than the noise (after 2x2 binning, the observed profiles have a 
signal-to-noise ratio of about 120 in the continuum intensity). 
Figure[3]shows examples of profiles recorded in the center- and 
limb-side penumbra, along with the best-fit profiles resulting 
from the inversion. The very asymmetrical shapes induced by 
the photospheric Evershed flow are successfully reproduced, 
indicating that the assumption of linear LOS velocity stratifica- 
tions suffices to explain the observed Doppler shifts and line bi- 
sectors. This is true for all profiles showing linear bisectors. In 
the very outer penumbra, bisector kinks and reversals are com- 
mon (Paper II). We cannot reproduce these profiles so well. In 
these cases, it seems that more complex LOS velocity stratifi- 
cations than the one used here would apply. 

3. Results 

The filtergrams have been binned by a factor of 2 before ex- 
tracting the intensity profiles. This increases the pixel size to 



0718 x 0718, which is still sufficient to sample our seeing- 
limited spatial resolution of about 075. After binning, we are 
left with 37 000 individual profiles which are inverted indepen- 
dently. The inversion has been performed using a parallelized 
version of SIR in a 16-processor Linux Beowulf cluster. The 
inversion of the data cube is accomplished in about 2.5 hours. 

Figure [4] displays maps of the temperature and LOS ve- 
locity at two representative optical depths in the atmosphere 
(log T500 = and -2). The fine structure of the penumbra is 
apparent (both in temperature and velocity) in the deep photo- 
sphere, where hot penumbral filaments and flow filaments are 
seen to extend from the inner to the outer penumbral boundary. 
The fluctuations in temperature and LOS velocity are much re- 
duced at log T500 = -2; indeed, no fine-scale organization of 
the penumbra is obvious at this optical depth. Consequently, 
the physical mechanisms producing the filaments must operate 
preferentially in the deep layers. 

The right panels of Fig. |4] show maps of the microtur- 
bulence at logrsoo = and the macroturbulent velocity. 
According to our results, no microturbulence is needed in either 
the umbra or the inner penumbra. This is in contrast to the outer 
penumbra, where microturbulent velocities of 1 .8-2 km s' 1 are 
common. Thus, the smaller line widths observed in the inner 
penumbra (Fig. 5 of Paper I) can be interpreted as being due 
to a reduced microturbulence. Note that the map of microtur- 
bulence shows fine structure in the mid penumbra. Also impor- 
tant is the fact that microturbulence decreases with height, until 
zero values are reached in the high layers (see Sect. 15. 3> . Again, 
this is consistent with the idea that the mechanisms enhancing 
the line width operate preferentially in deep atmospheric lay- 
ers. The macroturbulence stays relatively constant across the 
penumbra, with values of 1-1.2 kms . Only in the very in- 
ner penumbra and the umbra is the macroturbulence reduced, 
perhaps as a consequence of their stronger magnetic fields. 

In the next sections, we describe in more detail the thermal 
and kinematic configuration of the penumbra as deduced from 
the inversion. 

4. Thermal structure 

First, let us examine the reliability of the inferred temperatures. 
Figure [5^ shows the temperature at T557.6 = cos6> = 0.92 re- 
turned by the inversion code versus the observed continuum 
intensities. The scatter of the points is small, indicating an ex- 
cellent correlation between the two quantities. The fact that we 
use an absolute normalization for the observed profiles is one 
of the main reasons why the inversion code is able to determine 
the temperature of the continuum-forming layers so accurately. 

The solid line in Fig. [5^ represents the temperatures ob- 
tained from the Eddington-Barbier approximation h(9) = 
S a(i~a - cos 8), where h(9) is the specific intensity observed 
at an heliocentric angle 6, Sa is the source function (assumed 
to be the Planck function), and ta is the continuum optical 
depth at the wavelength of the observation. As can be seen, 
the Eddington-Barbier approximation yields the same temper- 
atures as those deduced from the inversion. Only for small con- 
tinuum intensities are the inferred temperatures slightly larger 
than those indicated by the Eddington-Barbier relation, but 
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Fig. 4. Maps of physical quantities at logrsoo = (top) and logrsoo = -2 (bottom). Depicted from left to right are the temper- 
ature, the LOS velocity and the microturbulence. Negative LOS velocities indicate redshifts. The lower right panel displays the 
macroturbulent velocity, assumed to be height independent. Contours outline the inner and outer penumbral boundaries. 



overall the agreement is very satisfactory. This implies that the 
observed continuum intensities can be reliably used to estimate 
temperatures in deep photospheric layers. A similar conclusion 
has been reached by TTR and Rouppe van der Voort (2002). 

In the line wing, the one-to-one correspondence between 
specific intensity and temperature found before does not hold 
any longer. This is demonstrated in Fig.|5j3, where we plot the 
temperatures resulting from the inversion at logrsoo = —1.3 
versus the intensities observed at -5.2 pm from line center. We 
use logTsoo = - l-3 because this is the optical depth of the 
centroid of the contribution function for the line depression at 
AA = -5.2 pm, evaluated in the mean penumbral model atmo- 
sphere. It is immediately apparent from the figure that the same 
line wing intensity can be associated with a wide range of tem- 
peratures (up to about 500 K). Thus, line wing intensities can- 
not be employed to infer temperatures directly. The reason is 
that the intensities critically depend on the line asymmetries in- 



duced by the Evershed flow and on the line width variations due 
to microturbulent and macroturbulent velocities. All these fac- 
tors contribute to the large scatter of the points seen in Fig. \5]p. 
The problem has been recognized by Balasubramaniam (2002). 
In an attempt to cure it, Balasubramaniam introduced the con- 
cept of flowless maps, i.e., intensity maps constructed by re- 
moving the global Doppler shift of the lines emerging from 
each individual pixel. Unfortunately, flowless maps are still af- 
fected by the line asymmetries caused by the Evershed flow. 

Let us now investigate the radial variation of the tempera- 
ture at different optical depths. To this end, we have computed 
azimuthally averaged temperatures along elliptical paths cen- 
tered on the spot (cf. Fig. 6 of Paper I). The average tempera- 
tures are plotted in Fig.[6]for four optical depths. At all heights 
in the atmosphere, the temperature is observed to increase radi- 
ally from the sunspot center to the outer penumbral boundary. 
This increase is rather linear except for an obvious hump in the 
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Fig. 5. (a) Temperature at T557 6 = cos 6 = 0.92 inferred from 
the inversion vs continuum intensity (normalized to the con- 
tinuum of the average quiet sun profile). The solid line rep- 
resents the temperatures obtained from the Eddington-Barbier 
approximation using the observed continuum intensities, (b) 
Temperature at logr 5 oo = -1.3 vs the intensity observed at 
A/l = -5.2 pm from line center, normalized to the continuum 
of the average quiet sun profile. 



inner penumbra, where the temperature seems to be enhanced 
by several hundreds K. The amplitude of the hump decreases 
as one moves towards the upper photospheric layers. In deep 
layers (log T500 = 0), the hump is so pronounced that the ra- 
dial variation of the temperature is almost flat between 0.6 and 
0.9 penumbral radii (i.e., the mid penumbra and most of the 
outer penumbra). We also note that the hump is located increas- 
ingly closer to the umbra/penumbra boundary as the upper lay- 
ers are approached. These temperature enhancements could be 
due to the presence of hot penumbral tubes, as suggested by 
the numerical simulations of Schlichenmaier et al. (1998). The 
hot tubes would be more easily detected in the inner penum- 
bra because the plasma emerging from subphotospheric layers 
quickly cools off by radiation away from the tube's inner foot- 
point (Schlichenmaier et al. 1999). If the humps exhibited by 
the temperature curves of Fig. [6] are indeed produced by hot 
flux tubes, then it is clear that they must be located preferen- 
tially in the lower layers. 

The mean temperature stratification of the penumbra con- 
sidering all pixels is shown in Fig. For comparison, we 
also plot the mean temperature stratification obtained by TTR. 
The two curves have the same curvature at all heights, but are 
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Fig. 6. Radial variation of the azimuthally averaged tempera- 
tures at several optical depths. The shaded areas represent the 
standard deviation of the temperatures at each radial distance. 
The vertical line marks the inner penumbral boundary. 
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Fig. 7. (a) Average temperature stratification in the penumbra 
(solid). The penumbral model of TTR (shifted by 200 K) is also 
shown for comparison (dashed). The shaded areas represent the 
rms fluctuations of the temperature, (b) Average temperature 
stratification in the center-side and limb-side penumbra (solid 
and dashed, respectively). 

shifted by 200 K (our temperatures being cooler than those of 
TTR). We do not deem that the offset of 200 K is significant, 
as different spots may certainly show different temperatures. In 
addition, the two data sets may be affected by different levels 
of stray light contamination. 

In Fig. 05 we plot the average temperature stratifications 
for the center-side and limb-side penumbra. Again, the curva- 
ture of the two stratifications is very similar, but the limb-side 
penumbra is observed to be cooler than the center-side penum- 
bra at all heights. The difference is larger in the upper layers (up 
to 150 K). A similar asymmetry has been reported by Rouppe 
van der Voort (2002) from an analysis of Can K observations. 
Rouppe van der Voort suggests that the larger average temper- 
ature of the center-side penumbra is produced by an excess of 
bright structures in the outer center-side penumbra. Our results 
support this idea. In Fig. [8] histograms of the temperature at 
log T500 = are presented for the inner and outer parts of the 
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(solid and dashed lines, respectively). 



center and limb side penumbra. The inner penumbra does not 
show significant differences between the two sides of the spot. 
Pixels in the outer penumbra, by contrast, are generally hotter 
on the center side. 

The fact that both sides of the penumbra show different 
temperatures might be related to the different viewing angle. 
Rouppe van der Voort (2002) suggests that the higher tem- 
peratures of the outer center-side penumbra can be explained 
by isotherms being tilted downwards away from the umbra. 
The required tilt is opposite to what one would expect from 
the Wilson depression (inclined upwards away from the um- 
bra, e.g., Mathew et al. 2004), but it may be real: it has been 
demonstrated that penumbral flux tubes dive down below the 
solar surface in the mid penumbra and beyond (Westendorp 
Plaza et al. 1997, 2001; Bellot Rubio et al. 2004; Borrero et al. 
2004). That is, the field lines are slightly inclined downwards in 
the outer penumbra. Perhaps this configuration of the magnetic 
field is capable of determining the inclination of the isotherms, 
although the exact mechanism has not been identified yet. 

As pointed out by Rouppe van der Voort (2002), other au- 
thors also find that the center-side penumbra is generally hotter 
than the limb-side penumbra (e.g., TTR; Westendorp Plaza et 
al. 2001). To the best of our knowledge, only one discrepant 
analysis exists. Schmidt & Fritz (2004) have studied the az- 
imuthal variation of the continuum intensity in a number of 
spots observed at different heliocentric angles. These authors 



log t 500 = 




100 



200 



300 




100 



200 300 
Azimuthal angle [deg] 



Fig. 9. Temperature fluctuations at different optical depths 
along an azimuthal path at r/R = 0.8. The vertical dotted lines 
mark a number of selected positions showing enhanced temper- 
atures. Azimuthal angles are measured counterclockwise, with 
90° pointing to the limb and 270° to disk center. 



found larger intensities in the limb-side penumbra for spots 
close to the solar limb, i.e., the opposite behavior. Clearly, more 
observations are required to settle the issue. 

Taking advantage of the high spatial resolution of our 
dataset, we conclude this section with a discussion of the 
temperature fluctuations observed at different optical depths. 
Figure [5] displays such fluctuations along an azimuthal path 
crossing the penumbra at a normalized radial distance of 0.8, 
for logTsoo = 0, -1, and -2. The fluctuations become smaller 
and smaller as the upper layers are approached. At log T500 = 0, 
the temperature fluctuates by some 200 K, while at log T500 = 
-2 the fluctuations are smaller than 50-100 K. The vertical 
lines mark a few positions where the temperature is enhanced 
locally. Clearly, hot structures in the deep layers remain hot 
in the upper layers, but with lower temperature contrasts. This 
is in excellent agreement with the findings of Rouppe van der 
Voort (2002) and Bello Gonzalez et al. (2005). Our results sug- 
gest that the structures giving rise to the bright penumbral fila- 
ments seen in white-light images (probably flux tubes) are lo- 
cated mainly in the lower photosphere. The small temperature 
fluctuations observed in high layers do not necessarily imply 
that flux tubes are present there: if the deep-lying tubes are hot, 
it is very likely that they can heat the layers above them. Thus, 
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Fig. 10. Comparison of LOS velocities resulting from the in- 
version and bisector velocities at intensity levels of 60%-80% 
(left) and 2%-20% (right). The solid lines represent a one-to- 
one correspondence. 



small temperature fluctuations in high layers could be produced 
by tubes located deeper down. 



5. Kinematic structure 

We start this section by examining how well the bisector veloci- 
ties of Papers I and II compare with the velocities inferred from 
the inversion. Such a comparison is important to strengthen the 
conclusions drawn in Paper II from the analysis of the observed 
bisectors. Figure [TUI shows scatter plots of the LOS velocities 
deduced from the two different methods. As can be seen, the 
agreement is very satisfactory in both deep and high layers. 
Bisector velocities near the continuum (intensity levels 60%- 
80%) are well correlated with the LOS velocities determined 
at logrgno = -0.6, with an rms difference of only 55 ms -1 . 
Bisector velocities near the line core (intensity levels 2%-20%) 
are also well correlated with the LOS velocities returned by the 
inversion at logrson ~ -1.8, the rms difference being about 
70 m s _1 . We have chosen the optical depths that yield the best 
agreement, but it should be remarked that they coincide almost 
exactly with the centroids of the RFs to velocity perturbations, 
evaluated in the average penumbral model atmosphere. The 
centroids are located at logrsno 0.8 and ~ -1.6 for inten- 
sity levels of 60%-80% and 2%-20%, respectively. The agree- 
ment is not surprising: Ruiz Cobo & del Toro Iniesta (1994) 
and Sanchez Almeida et al. (1996) have shown that physical 
parameters measured from spectral lines using simple tech- 
niques (e.g., Doppler shifts from line bisectors) are an average 
of the actual stratification of the parameter weighted by the cor- 
responding RF. If the stratification is linear with depth, then the 
height of formation of the measured parameter is the bary center 
of the RF (del Toro Iniesta 2003, Chapter 10). 

Figure ^3 demonstrates that Doppler velocities derived 
from bisectors near the line core should be ascribed to layers 
significantly deeper than the 'formation height' of the line core. 
This is because of the very large widths of the RFs of / to v L qs- 
Such widths imply that a wide range of layers contribute to the 
observed Doppler shift. As shown in Fig. 1, the line core reacts 
only little to mass motions in the upper photosphere, and so 
velocities there are difficult to detect unless the contribution of 
deeper layers is properly considered. 



5.1. Height variation of the LOS velocity 

In Fig.^2 we show the azimuthally averaged LOS velocity at 
four optical depths as a function of normalized radial distance. 
We consider separately the limb-side and center-side penum- 
bra (left and right upper panels, respectively). The curves de- 
picted in Fig. II llconfirm in a more quantitative way the impres- 
sion from Fig.|4]that the Evershed flow strongly decreases with 
height. A velocity reversal is observed around log T500 = -2. 
On the limb side, we infer redshifts in the deep photosphere 
and blueshifts in the upper layers. For most of the center-side 
penumbra, blueshifts in deep layers and redshifts in high layers 
are retrieved (an exception is the very outer part of the center- 
side penumbra). Interestingly, the LOS velocity stratifications 
displayed in Fig. 4 of TTR also feature sign reversals around 
log T500 = -2. A possible explanation for the sign reversal is 
offered in Sect. l6.2l 

The lower panels of Fig.^2 snow the average velocity strat- 
ifications at two radial distances: r = 0.6R and r = Q.91R. On 
the limb side, the gradient of vlos with logrson is significantly 
larger in the outer penumbra. On the center side, the situation is 
reversed, with slightly larger gradients in the inner penumbra. 
These gradients reproduce the slopes of the observed bisectors 
in the different parts of the penumbra (Fig. 1 of Paper II). Note 
the occurrence of negative velocities (blueshifts) in the upper 
layers except in the inner center-side penumbra, where posi- 
tive velocities are retrieved above logrsoo = -2. The inferred 
velocity stratifications, together with the various broadening 
and smearing mechanisms, explain the line-core blueshifts ob- 
served almost everywhere in the penumbra, as well as the small 
line-core redshifts detected in the inner center-side penumbra 
(e.g., bisectors #33 and #34 in Paper II). At this point it is im- 
portant to remark that the line-core blueshifts cannot be pro- 
duced by the inverse Evershed flow. According to Georgakilas 
et al. (2003, his Fig. 6), the speed and angle with respect to the 
local vertical of the inverse Evershed flow are ~ 0.5 km s _1 and 
100° in the inner penumbra, at the height where the intensity 
observed at 0.5 A from the center of the Ha line is formed. 
The flow speed and the flow angle near the outer penumbral 
boundary are ~ 4.5kms _1 and 115°, respectively. With these 
flow configurations and the heliocentric angle of our observa- 
tions, the inverse Evershed effect would produce non-negligible 
Doppler shifts only in the outer center-side penumbra, where it 
would show up as a redshift, not as a blueshift. 

5.2. Flow geometry at different heights 

We have computed the flow speed and flow angle at different 
optical depths from the azimuthal variation of the LOS veloci- 
ties returned by the inversion code. The calculations have been 
performed as explained in Sect. 5.4 of Paper I. If one assumes 
that the flow field is axially symmetric, the azimuthally aver- 
aged LOS velocity reflects the vertical flow component, and 
the amplitude of the azimuthal variation of the LOS velocity re- 
flects the horizontal flow component. Using these components, 
the mean flow inclination and the absolute flow velocity can be 
deduced. The results of the analysis are presented in Fig.ll2las 
a function of normalized radial distance. 
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Fig. 11. Top: Radial variation of the azimuthally averaged LOS velocity in the limb-side (left) and center-side penumbra (right), 
for four different optical depths. The shaded areas represent the rms variation of the LOS velocities at each radial distance. 
Positive LOS velocities indicate motions away from the observer (i.e., redshifts). Bottom: Average LOS velocity stratifications 
in the limb-side penumbra (left) and center-side penumbra (right) at r — 0.6R and r = 0.97 R. 



As mentioned in Sect. 12.31 our simple model atmosphere 
does not account for the possibility of unresolved structure in 
the resolution element. Thus, the flow speeds resulting from the 
calculations are probably smaller than the real ones 1 . To first 
order, however, the flow angles are not influenced by this effect, 
because they are derived from the ratio of the horizontal and 
vertical components of the flow, and both are affected equally 
by the lack of spatial resolution. 

From Fig.[21we note the following features: 

1 . The flow inclination (with respect to the local vertical) in- 
creases with optical depth. Downflows are present in the 
mid and outer penumbra, but only in deep layers. In high 
layers, the flow is always inclined upward. 

2. The flow velocity decreases with decreasing optical depth. 

3. The radial position of maximum flow velocity and maxi- 
mum flow inclination 'migrate' outwards as we go to higher 
atmospheric layers. 



Inclination angle 



Flow velocity [km/s] 



1 Following Bellot Rubio (2004), we may assume for simplicity that 
at a given radial distance the flow channels occupy a fraction < a < 
1 of the resolution element, the remaining (1 - a) being plasma at 
rest. As a rough estimate, the LOS velocity observed in the resolution 
element would be vlos ~ ffvixis, with vlos the intrinsic LOS velocity 
of the flow channel. 
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Fig. 12. Flow geometry at different optical depths as deduced 
from the azimuthal variation of the observed LOS velocity. 
Left: Flow inclination (with respect to the local vertical). Right: 
Absolute magnitude of the flow. 



In summary, the azimuthally averaged flow velocity in- 
creases with depth and the azimuthally averaged flow inclina- 
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tion show downflows only in the deep layers. This strength- 
ens the conclusions of Paper II, where we inferred a similar 
flow configuration from the interpretation of individual bisec- 
tor shapes. It is important to remark that calculations using the 
bisector velocities derived in Paper II lead to the very same 
flow angles and flow speeds (cf. Fig. 6 in Bellot Rubio 2004), 
as could have been expected from Fig. 1101 

The results presented in Fig.^]are in good agreement with 
the flow inclinations determined by Bellot Rubio et al. (2003, 
2004) from a two-component inversion of infrared Stokes pro- 
files. The radial variation of the flow angle described by Bellot 
Rubio et al. corresponds roughly to our curve for log T500 = 0, 
because the Fei 1565 nm lines they use are sensitive to veloc- 
ities in very deep and narrow photospheric layers. Our results 
are also consistent with those of Sanchez Cuberes et al. (2005), 
who inverted spectropolarimetric measurements of a sunspot at 
disk center in terms of one-component model atmospheres with 
gradients of velocity. These authors found that downflows oc- 
cur predominantly in deep layers (log T500 < -1), disappearing 
higher up in the atmosphere. The present analysis, however, 
does not confirm the small downflows retrieved by Sanchez 
Cuberes et al. (2005) in the inner penumbra. 




Normalized radial distance 

Fig. 13. Radial variation of the azimuthally averaged micro- 
turbulence at four optical depths. The shaded areas represent 
the standard deviation of the microturbulent velocities found at 
each radial distance. The vertical line marks the inner penum- 
bral boundary. 



5.3. Microturbulence 

For each pixel the inversion returns the height variation of the 
microturbulence, assumed to be linear with log T500. Figure [T31 
shows the azimuthally averaged microturbulence at four optical 
depths versus normalized radial distance. From this figure it is 
apparent that: 

1 . In the inner penumbra, up to about 0.6 penumbral radii, the 
microturbulence is zero at all heights. 

2. In the mid and outer penumbra, the microturbulence is dif- 
ferent from zero and increases with optical depth and radial 
distance. 

3. The gradient with depth of the microturbulence increases 
with increasing radial distance. 

In Paper I (Fig. 7) we determined the radial dependence of 
the azimuthally averaged equivalent width (EW) and the full 
width at half maximum (FWHM) of the intensity profiles ob- 
served in NOAA 10019. For both parameters we found minima 
in the inner penumbra (r ~ 0.65/?) and maxima in the outer 
penumbra (r ~ Q.95R). According to the present results, this 
radial increase of the EW and the FWHM is due to a microtur- 
bulence that increases radially (item 2 above). 

What is the origin of the microturbulent velocities re- 
turned by the inversion code? One can speculate that stronger 
Evershed flows result in increased shear at the boundary layer 
separating the flow channels from the surrounding (static) mag- 
netic atmosphere. Strong shears may produce small-scale tur- 
bulence through the onset of instabilities. This would explain 
why we detect enhanced microturbulence in the mid and outer 
penumbra, where the flow speeds are larger. 

However, we also found in Paper I that both the EW and 
the FWHM tend to be slightly enhanced in the regions of max- 
imum Doppler shifts (i.e., along the line connecting the disk 
center and the sunspot center). A similar tendency is observed 




LOS velocity at log x 500 = -1 .5 [km s' 1 ] 

Fig. 14. LOS velocities determined from the inversion at 
logrsoo = -1-5 versus the microturbulence at the same opti- 
cal depth. The vertical dotted lines indicate LOS velocities of 
±0.4 km s -1 . Larger velocities are always associated with en- 
hanced microturbulence. The same correlation is seen in deeper 
layers, but less clearly. 



for the microturbulence. Figure [2] shows the microturbulence 
determined from the inversion at logrsoo = -1-5 versus the 
LOS velocity at the same optical depth. From this figure, it 
is clear that large LOS velocities are always associated with 
large microturbulences. The findings of Paper I led us to con- 
jecture that the increase in the EW and the FWHM associated 
with large velocities may be due to the presence of strongly 
Doppler-shifted (but spectrally unresolved) line satellites. This 
suggests that part of the microturbulence we infer may not rep- 
resent real small-scale turbulence, but a convenient means to 
account for the effects of line satellites. 
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6. Discussion 

6.1. Origin of enhanced line-wing intensities in the 
inner penumbra 

From an analysis of the Stokes profiles of the Fei 1564.8 nm 
line, Bellot Rubio (2003) found that the intensity emerging 
from the inner penumbra is higher than that emerging from the 
outer penumbra at wavelengths close to, but not exactly at, the 
line center. That is, sunspots exhibit a bright ring in the in- 
ner penumbra when observed in the wing of spectral lines. The 
bright ring can be interpreted as a temperature enhancement in 
the mid-photospheric layers (as it is detected in the wings of 
the line but not in the continuum or line core). 

Our high spatial resolution observations of NOAA 10019 
confirm this behavior (Fig. 8 of Paper I) and at the same time 
add new information on it: the intensity enhancement occurs 
in the form of radial, narrow fibrils very similar to the bright 
penumbral filaments seen in continuum images, but with much 
higher contrast (Fig. 3 of Paper I). Moreover, the enhancement 
is observed preferentially in the inner center-side penumbra 
blueward of the line core, and in the inner limb-side penum- 
bra redward of the line core. The temperature stratifications 
of Fig. |6] indicate that the inner penumbra is always cooler 
than the outer penumbra at all heights in the atmosphere, so 
the enhanced brightness cannot be due to a temperature effect. 
Heating of the mid photosphere was the mechanism originally 
proposed by Bellot Rubio (2003). It is also the scenario favored 
by Sanchez Cuberes et al. (2005). 

Here we argue that the bright ring is actually produced by 
(a) the small Doppler shifts and line asymmetries observed 
in the inner penumbra, together with (b) the increase of the 
line width toward the outer sunspot boundary. The left panel 
of Fig. 1151 shows an intensity map of the spot at AA = -6.4 
pm from line center (filtergram #42 in Fig. 3 of Paper I). The 
enhanced brightness of the inner center-side penumbra is ap- 
parent in this image. The intensity profiles emerging from the 
two spatial positions marked with crosses are displayed in the 
right panel of Fig. [H] As can be seen, these profiles possess 
very different properties: in the outer penumbra (solid line) 
they are broader and have a stronger asymmetry than in the 
inner penumbra (dashed line). The bright ring results from a 
comparison of the intensities of pixels belonging to the inner 
and outer penumbra at fixed wavelengths (dotted vertical lines). 
Consider, for example, wavelengths in the blue wing (leftmost 
vertical line). At this wavelength, the inner penumbra is clearly 
brighter than the outer penumbra. In the red wing (rightmost 
vertical line), the inner penumbra would still be brighter, but 
with a much reduced contrast. This is consistent with the ob- 
servations described in Paper I. Similar arguments apply to pro- 
files emerging from the limb-side penumbra: because of the 
smaller line widths and asymmetries in the inner penumbra, the 
intensity at fixed wavelengths is enhanced relative to the outer 
penumbra, with higher contrasts redward of the line core. 

Thus, the bright ring can be explained satisfactorily in 
terms of the radial variation of the microturbulence and the 
Evershed flow. Microturbulence is an important parameter: if it 
were the same everywhere across the spot, redward of the line 



i 

_Z7 



Center-side profiles 




-30 -20 -10 10 20 30 
AX [pm] 

Fig. 15. Left: Intensity map of the spot at A/l = -6.4 pm from 
line center. Two pixels belonging to the inner and outer center- 
side penumbra are marked with crosses. Right: Best-fit inten- 
sity profiles returned by the inversion at the spatial positions 
marked with crosses in the left panel. Solid: outer penumbra. 
Dashed: inner penumbra. The vertical dotted lines mark fixed 
wavelength positions in the blue wing, the line core, and the 
red wing (from left to right). The blue and red wing positions 
correspond to filtergrams #42 and #58 of Fig. 3 in Paper I. 
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Fig. 16. Maps of temperature at logrsoo = (left) and 
logTsoo = -1 (right) resulting from an inversion of the ob- 
served profiles in which only constant velocities along the LOS 
are permitted. 



core the inner penumbra would appear darker than the outer 
penumbra on the center side, contradicting the observations. 
To further support the idea that the bright ring is not due to 
temperature enhancements, we have carried out the same in- 
versions by forcing the LOS velocity to be constant with height 
(i.e., we do not allow for line asymmetries). In addition, we set 
the microturbulence to zero. The thermal stratification resulting 
from these inversions is shown in Fig. [HO at logTsoo = and 
-1. While the temperature in the deep layers is essentially un- 
changed, the map of temperature at log T500 = — 1 shows a nice 
bright ring in the inner penumbra. From this we conclude that 
the only way to explain the enhanced line-wing intensities of 
the inner penumbra when no line asymmetries are considered 
is by an increase of the temperature in the mid photosphere. 

Using flowless maps, Balasubramaniam (2002) detected 
bright structures in the inner penumbra at wavelengths close to 
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line center (his Fig. 8). Transforming observed intensities into 
radiation temperatures, he concluded that these structures are 
hotter than their surroundings. In our opinion, however, they 
more likely trace regions of reduced Evershed flows and small- 
scale turbulence. 



6.2. Interpretation of the LOS velocity stratifications 

In Paper II we demonstrated that the uncombed model pro- 
posed by Solanki & Montavon (1993) is able to explain the 
bisector shapes observed in the penumbra of NOAA 10019. In 
an uncombed penumbra, the Evershed flow occurs along flux 
tubes whose magnetic fields are more horizontal than that of 
their surroundings. Here we use this model to interpret the LOS 
velocity stratifications displayed in Fig. II II 

As pointed out in Sect. 12.31 it is the magnitude and slope 
of the inferred velocity stratification which inform us about the 
flow geometry (flow speed and flow angle) and height of the 
Evershed channels. Figure 5 of Paper II shows that, in the very 
outer penumbra, two flow channels stacked along the LOS and 
with different inclination angles are necessary to understand 
the kinks and reversals of the observed bisectors. In the inner 
and mid penumbra the bisectors are rather linear, and can be 
explained in terms of a single, deep-lying flow channel. 

Based on these results, in the inner penumbra (r = 0.6/?) 
we consider a flow channel located in deep layers, between 
logTgoo = -0.7 and 0. The flow has a speed of 8 km s _1 and 
an inclination of 87° with respect to the vertical. In the outer 
penumbra (r = 0.97/?), we assume two flow channels located 
in deep and mid photospheric layers. The deep-lying channel 
extends from logT5oo = -0.2 to 0.5. It has a flow speed of 8 
km s _1 and an inclination of 135°. The upper channel extends 
from logTsoo = —1.8 to —1.1 and is assumed to have a flow 
speed of 7 km s _1 and an inclination of 85°. This flow config- 
uration is similar to the one deduced in Paper II, except for the 
fact that our higher-lying channel is slightly inclined upward 
instead of being horizontal. The LOS velocities resulting from 
the two flux-tube geometries at the heliocentric angle of our 
observations (6 = 23°) are represented by solid lines in Fig. 1171 
for the center and limb-side penumbra. 

To investigate whether or not these flow configurations can 
explain the LOS velocities retrieved from the inversion, we 
have computed synthetic profiles under the assumption that 
the resolution element contains two different atmospheres. The 
temperature and pressure stratifications of both atmospheres 
are taken to be the mean stratifications deduced from the in- 
version of NOAA 10019. The LOS velocity is set to zero in 
the first component, whereas for the second component we use 
the stratifications of Fig. In the inner penumbra we adopt 
a filling factor of 0.2, i.e., only 20% of the resolution element 
is occupied by flows. In the outer penumbra the filling factor 
is 0.5, as in Paper II. The resulting profiles are convolved with 
the instrumental profile of TESOS and a macroturbulent veloc- 
ity of 2 km s _1 (inner penumbra) or 1 km s _1 (outer penumbra). 
Finally, they are inverted under the same conditions as the real 
data (two nodes for velocity, temperature and microturbulence, 
and a height-independent macroturbulence). Not unexpectedly, 
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Fig. 17. LOS velocity stratifications due to flux tubes along the 
LOS (solid lines) that would result in linear velocity stratifica- 
tions (dashed lines) when two nodes are used to retrieve the ac- 
tual velocity stratification. See text for details on the flux-tube 
geometries adopted in the inner and outer penumbra. 



the simulated profiles show the same bisector shapes and bisec- 
tor variations across the penumbra as the real observations. 

The linear velocity stratifications retrieved from the inver- 
sion of the simulated profiles are represented by dashed lines 
in Fig.^] The first thing to notice is that the inferred stratifi- 
cations change sign in the mid photosphere (around log T500 = 
-2) both on the center and limb sides, even though there are 
no structures in high layers that could be responsible for such 
velocity reversals. Actually, the sign reversal is the result of ex- 
trapolating the jump(s) of velocity sensed by the LOS to the up- 
per layers, where the response of the line to plasma motions is 
small. This extrapolation effect was first described by Martinez 
Pillet (2000) in connection with inversions of Stokes profiles 
of visible lines that do not explicitly account for the existence 
of penumbral flux tubes. The same effect has been observed 
by Mathew et al. (2003) in their one-component inversions of 
near-infrared lines. 

The second thing to note is that the azimuthally averaged 
LOS velocities inferred from the inversion of NOAA 10019 
(Fig. II 1> are in excellent accord with the linear velocity strat- 
ifications resulting from the flux tube configurations displayed 
in Fig. Specific features that are well reproduced include 
the following: 

- On the limb side, the slope of the stratifications is larger in 
the outer penumbra as compared with the inner penumbra, 
whereas on the center side it is slightly larger in the inner 
penumbra. 

- The slope of the stratifications in the inner penumbra is 
larger on the center side. 
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Fig. 18. Continuum intensity fluctuations (solid) and LOS velocity at logTsoo = (dashed) along an azimuthal path crossing the 
inner center-side penumbra as indicated on the left image. Negative velocities represent blueshifts. The vertical dotted lines mark 
positions of local maxima of the continuum intensity. 



- The LOS velocities corresponding to the outer center-side 
penumbra are negative all along the atmosphere. 

From this we conclude that the presence of flux tubes in the 
penumbra can explain the observed LOS velocity stratifica- 
tions. In particular, the seemingly weird run with depth of the 
LOS velocity in the center-side penumbra at normalized ra- 
dial distances larger than 0.9 (with zero or small redshifts at 
log T500 = and blueshifts increasing with height, see Fig. II 1> 
can be understood in terms of two flux tubes stacked along the 
LOS. The tube lying deeper contributes with positive veloc- 
ity (i.e., a redshift). This redshift is necessary to explain the 
bisector reversals described in Paper II. In order to produce a 
redshift in the center-side penumbra, the deep-lying flux tube 
must be inclined downward with respect to the horizontal by 
more than 23° (the heliocentric angle of the spot), otherwise it 
would show negative velocities and the slope of the resulting 
LOS velocity stratification would be opposite to what we ob- 
serve. Flux tubes returning to the solar surface have previously 
been inferred by, among others, Westendorp Plaza et al. (1997, 
2001), Schlichenmaier & Schmidt (2000), Bellot Rubio et al. 
(2003, 2004) and Borrero et al. (2004, 2005). 

The tube located higher in the atmosphere may be horizon- 
tal (as in Paper II) or slightly inclined upwards (less than 23° 
from the horizontal). The later configuration would still pro- 
duce redshifts on the limb side and blueshifts on the center side, 
exactly what is required to explain the velocity stratifications 
observed in the penumbra. Interestingly, such high-lying tubes 
do not seem to be necessary in the inner penumbra, only in the 
mid and outer penumbra. It is tempting to associate them with 
those described by Westendorp Plaza et al. (2001) and Bellot 
Rubio et al. (2004). These authors found indications of a new 
family of penumbral tubes in the mid penumbra and beyond. 
Spectropolarimetric measurements suggest that such tubes are 
more vertically oriented than the deeper ones, carrying a frac- 
tion of the Evershed flow to the upper photosphere. 

6.3. Do we see individual flow channels? 

Continuum intensity fluctuations and LOS velocities are very 
well correlated in the inner penumbra. This is shown in Fig. 1181 
where we plot normalized continuum intensity variations and 
LOS velocities at log t 5 qo - along an azimuthal path crossing 



8 distinct continuum filaments of the inner center-side penum- 
bra (we do not consider the limb-side part of the azimuthal path 
because no clear filaments are observed on the limb side). The 
curves depicted in Fig.^Dreveal that larger continuum intensi- 
ties are associated with larger LOS velocities, i.e., the Evershed 
flow is more intense in the brighter structures, at least in the in- 
ner center-side penumbra. This agrees with the conclusions of 
Bello Gonzalez et al. (2005) and Schlichenmaier et al. (2005). 

However, the most important fact demonstrated by Fig. 1181 
is that also the dark structures show large LOS velocities. Stray 
light from the bright filaments may change the line profile in 
the dark structures, but only to some extent (cf. Sect. 3.2 in 
Paper I). Unless the stray light contamination is unreasonably 
large, the velocities found in the dark filaments must be real. 
This means that we see the Evershed flow everywhere: even at 
a resolution of (J.'5, there is almost no position in the penumbra 
where the velocity drops to zero (except where it must be zero 
for geometrical reasons). In our opinion, this well established 
fact (see also Hirzberger & Kneer 2001 and Rouppe van der 
Voort 2002) has not been emphasized adequately in the past. 

Since Evershed flows are ubiquitous, it is tempting to think 
that the bright and dark structures observed in continuum im- 
ages are individual flow channels (i.e., isolated flux tubes). But 
this interpretation is problematic. For example, one would have 
to explain why the flow channels possess different tempera- 
tures and velocities at the same radial distance and, more im- 
portantly, why there exists such a regular azimuthal pattern of 
temperature enhancements/coolings and larger/smaller plasma 
flows. An alternative explanation is that the filaments seen in 
continuum images are not individual flow channels, but a col- 
lection of them. If the filling factor of the tubes changes regu- 
larly in the azimuthal direction, i.e., if the number of flow chan- 
nels present in the resolution element changes azimuthally, then 
one would observe larger LOS velocities and higher tempera- 
tures at those positions where there are more tubes (provided, 
of course, that the tubes are hot in the inner penumbra, as sug- 
gested by the simulations of Schlichenmaier et al. 1998 and 
the uncombed inversions of Borrero et al. 2005). The impor- 
tant idea here is that the intrinsic flow speed and temperature 
of the flux tubes would not change azimuthally at a fixed ra- 
dial distance from sunspot's center, only the filling factor would 
change. We believe this is a natural explanation for the correla- 
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tion displayed in Fig.^] Independent support for this idea has 
been gathered from the analysis of full Stokes profiles of near- 
infrared lines (Bellot Rubio 2004; Bellot Rubio et al. 2004). 

Our results strongly suggest that the filaments observed in 
Doppler maps at a resolution of about Of.' 5 are still unresolved. 
That is, they are not individual flow channels. For this reason, 
in Paper II we referred to them as flow filaments, to clearly 
distinguish them from the actual, smaller flow channels. 

7. Summary 

Using high spectral and spatial resolution intensity profiles of 
the non-magnetic Fei 557.6 nm line, we have studied the ther- 
mal and kinematic properties of a sunspot at 23° from disk cen- 
ter. The profiles, obtained with TESOS, have been inverted in 
terms of a one-component model atmosphere with gradients of 
the physical quantities. 

The temperature, LOS velocity, and microturbulent veloc- 
ity maps resulting from the inversion show that the fine struc- 
ture of the penumbra is less apparent in high layers than in 
deep layers. This suggests that the physical mechanism(s) pro- 
ducing the penumbral filaments operate mainly in the lower 
photosphere. 

Our analysis confirms the existence of a thermal asymmetry 
between the center and limb side penumbra, the former being 
hotter by 100-150 K on average. Such an asymmetry seems to 
be related to the presence of a larger number of hot structures 
in the outer center-side penumbra, but its nature is still unclear. 

We have also investigated the origin of the enhanced line- 
wing intensities exhibited by sunspots in the inner penumbra. 
Our results indicate that the bright ring is produced by the 
smaller line widths and Doppler shifts observed in the inner 
penumbra as compared with the outer penumbra. Thus, it does 
not reflect temperature enhancements in the mid photosphere. 

The inferred LOS velocities have been used to determine 
the geometry of the Evershed flow assuming axial symmetry 
of the flow field. We find that both the flow speed and flow an- 
gle increase with optical depth and radial distance. Downflows 
are detected in the mid and outer penumbra, but only in deep 
layers (logTsoo ^ -1.4). By means of simulated profiles, we 
have shown that the velocity stratifications retrieved from the 
inversion are consistent with the existence of penumbral flux 
tubes channeling the Evershed flow. In the inner penumbra, 
a deep-lying tube slightly inclined upwards is able to explain 
the observations. In the outer penumbra, two flux tubes stacked 
along the LOS and with different inclination angles are neces- 
sary. The deep-lying tube is inclined downwards, whereas the 
higher-lying tube carries upflows. We suggest that the later may 
belong to a different family of tubes, perhaps the ones first de- 
scribed by Westendorp Plaza et al. (2001). 

Finally, we have demonstrated that there is an excellent 
correlation between continuum intensities and LOS velocities 
in the inner center-side penumbra, with larger Evershed flows 
associated with brighter continuum structures. Very remark- 
ably, however, dark structures also exhibit significant Evershed 
flows. This led us to propose that the bright and dark filaments 
seen in continuum images at a resolution of 075 are not indi- 
vidual flow channels, but a collection of them. 



Our analysis shows that spectroscopy with 2D filter instru- 
ments provides a wealth of information on the physical proper- 
ties of sunspots. Obviously, the next step is to carry out 2D full 
spectropolarimetric observations near the diffraction limit. This 
will soon be possible with instruments like the Interferometric 
BIdimensional Spectrometer (IBIS) and the KIS/IAA Visible 
Imaging Polarimeter (VIP) attached to TESOS. 
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